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Abstract 

The nonlocality of the microscopic nucleon-nucleus optical potential is commonly localized by 
the Brieva-Rook approximation. The validity of the localization is tested for the proton+ 90 Zr 
scattering at the incident energies from 65 MeV to 800 MeV. The localization is valid in the wide 
incident-energy range. 
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I. INTRODUCTION 



Microscopic understanding of nucleon-nucleus (NA) elastic scattering is a long-standing 
fundamental subject in the nuclear reaction theory. This is nothing but to solve the many- 
body scattering problem. The many-body collision, however, can be approximately described 
as a scattering between two bodies interacting via a complex mean-field (optical) potential. 
This optical potential is an important ingredient in theoretical calculations of cross sections 
of elastic and inelastic scattering, charge exchange and transfer reactions, and so on. This 
means that a good global optical model is a powerful tool for predicting observables of NA 
scattering for which no measurements exist, e.g., scattering of unstable nuclei from proton 
target. 

A reasonable way of getting the optical potential is to calculate the NA folding potential 
with the nucleon-nucleon (NN) ^-matrix interaction . The interaction is first 

evaluated in infinite nuclear matter and then folded into target (A) density by using the 
local-density approximation. The ^-matrix interaction thus obtained is a complex nonlocal 
potential depending on the incident energy (E) of nucleon (N) and the nuclear-matter den- 
sity (p). This microscopic optical potential is successful in reproducing data of NA elastic 
scattering 4| in the wide range of 40 < E < 250 MeV from light to heavy targets. Above the 
pion production threshold, resonance and meson production effects are evident. Recently, a 
bare NN interaction was extrapolated to reproduce the NN scattering data to 2.5 GeV by 
adding a complex potential phenomenologically [3] , and the g- matrix interaction constructed 
from the complex NN interaction was also successful in reproducing the NA scattering at 
40 < E < 800 MeV 

In many applications, use of a nonlocal NA potential is impractical. For example, in 
8 B+A scattering the projectile easily breaks up into 7 Be and p. This projectile breakup 
processes are described by solving the scattering problem of three-body system 7 Be+p+A. If 
all potentials are local in the system, thisproblem can be solved by the method of continuum- 
discretized coupled channels (CDCC) [8|, 19]. If the potential between p and A and/or the 
potential between 7 Be and A is nonlocal, this is not easy. For such cases, use of an equivalent 
local potential is quite practical, if it is accurate. Brieva and Rook (BR) proposed an 
approximate form of the equivalent local potential [l|. This is commonly used in many 
applications; for example see Refs. {fj, [lo| and references therein. However, the validity of 
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the approximate form is not shown yet. 

In this paper, we show the validity of the BR localization over the wide range of 65 < 
E < 800 MeV, comparing the scattering solution of the non-local NA potential with that 
of the BR-type local potential. As a typical case, we consider the p+ 90 Zr scattering. The 
BR localization is composed of three approximations. We show that one of the three is 
redundant, and test the remaining two separately. 

In Sec. [Til the method of solving the Schrodinger equation with the non-local NA potential 
and the way of getting the ground-state wave function of target nucleus are presented. In 
Sec. IHIl the BR localization is recapitulated. In Sec. IIVI the validity of the BR localization 
and the related topics are argued. Section [V] is devoted to summary. 



II. FORMULATION 



In the g-matrix approach 



flfl 



6J, the microscopic NA optical potential is con- 
structed by folding the ^-matrix interaction g with the ground-state density of target 
nucleus (A), where S (T) is the spin (isospin) of the N+N system. In this procedure, the 
antisymmetrization between an incident nucleon and target nucleons in A is taken care of 
by using g ST which is properly antisymmetric with respect to the exchange of the colliding 



nucleons, since the prescription is shown to be a good approximation ll|, ll2j. In the ap- 
proximation, the folding potential is expressed by the sum of a local direct term £/ DR (R) 
and a ncnoca, exchange tenn ,-(R, r ) Q. Hence, the eiastic scattenn g can be described 
by solving the Schrodinger equation 

h 2 



--V 2 R + U^(R) + V c (R)5»_\ /2 



E 



Xk,^ (R) = J U EX (R, r) X K,, 1 (r) dr (1) 



for the relative wave function Xk,^ (R) ; where R stands for the coordinate of incident 
nucleon (N) from the center-of-mass of A, r is the coordinate of nucleon in A from the 
center-of-mass of A, V c (R) is the Coulomb potential, KK (E) is an incident momentum 
(energy), and V\ = 1/2 for neutron scattering and —1/2 for proton one. We assume that the 
target nucleus is much heavier than N. The relativistic kinematics is taken by defining the 
reduced mass jj as jj, = A/ m N + (p/c) 2 with p and mN the momentum and rest mass of N. 

In this paper, we consider only the central part of the microscopic optical potential, since 
it is a main component of the folding potential. We also assume that the ground state of 
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A is described by a single determinant of single-nucleon wave functions ^>u 2 ;nijj z ( r ,0> eacn 
classified with the z-component v<i of isospin, the principal quantum number n, the angular 
momentum I, and the total angular momentum j and its ^-component j z , where £ is the 
internal coordinate of the spin wave function 771/2 of nucleon. In this case, [7 DR (R) and 
U EX (R, r) are given by J U S, Q 

f/ DR (R) = W p V2 (r)g™(s;p U2 )dr, (2) 

f/ EX (R, r) = P«*( R > r )^f(«; P,J (3) 

with 

P^O) = Yl j V*u 2] nljjA r ^)Vu 2 ;nljjA r ^) d ^ ( 4 ) 
nljjz 

p V2 (R,r) = Y J V^nljj, ( r > Vmnljj, (R-, df, (5) 

<£ r =±i(*;p, 2 ) = ^WJ + S^WJ}^ (6) 
flgLofa P, 2 ) = ^ {^ 01 (a; P, 2 ) + 3s 10 (a; p„ 2 ) + (a; p, 2 ) + 3^ u (a; p„ 2 )} (7) 
^t z X = ± i P, 2 ) = ~ {V 1 (a; P, 2 ) + 3<? u (a; p, 2 )} (8) 

<?|f =0 (s;p, 2 ) = i{-^M S ;p, 2 )-3^ 10 ( S ;p, 2 )+^ 00 ( S ;p, 2 ) + 3^ 11 ( S ;p, 2 )}^ + ". ( 9 ) 
where s = r — R. The g matrix interaction g ST is a function of E and the single-particle 
density Pv 2 { r g)i where r g = \r g \ for the location r g at which the effective interaction works. As 
for the g matrix interaction, we take a sophisticated version of the Melbourne interaction 4| 
that is constructed from the Bonn-B NN potential [14| and includes a modification due to the 
pion-production effect Since the interaction has a finite range, r g can not be determined 
uniquely. Possible choices are (i) r g = r, (ii) r g = R and (iii) r g = r m = |r + R|/2. This 
ambiguity is referred to as the r 9 -ambiguity in this paper. The r 9 -ambiguity is small, as 
shown later in Fig. [HJ We then take choice (i). We will return to this point below. 
Expanding g® R into a series of multipoles, 



A 

one can get a simple form of 



<£ R («; p) = E 47r ^x~^ R A (r, p) (r) ® r A (? )] oo , (10) 



U DR (R) = 4vr W P» ( r ) S (r, R; P» 2 ) r 2 dr. (11) 

u 2 ,T z J 
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In this form, U BR (R) is a function of R = |R|. 

The scattering wave function XK,n( r ) is expanded into partial waves Xkv u l'{t) 



Air 



Xk, 



(12) 



L'M' 



Inserting Eq. (|T2l into Eq. (CD), multiplying the equation by Y^ M (R) from the left and 
integrating it over the solid angle R, one can get an equation for XKu lt L as 



E 



S2 



/2 



^ (L0A0|/0) 2 y 0: 2 . ny (r) sg A (r, i2; p„ 2 (r fl )) cf>„ 2 . nlj (R) X Ku u l (r) dr, (13) 

where gf x has been expanded into multipoles gf x A just as in Eq. (ITU]) , and (p U ;nij(r) is the 
radial part of the single-nucleon wave function l -Pv,nijjA r i £)■ m ^ ne derivation of Eq. ( {TBI . 

came out in 17 summation, so that Xkui,l(R) on the left hand side has the same angular 
momentum L as XK^^ij") on the right hand side. This is a consequence of the fact that 
Eq. ([TJ is rotational invariant in the coordinate space. In the present paper, Xkvi,l{R) is 
obtained by solving Eq. (|T3|) iteratively. In Eq. (fTB"]) . each multipole gf z x A depends on r g 
through p U2 , but it includes no information on an angle between vectors r and R. Hence, 
we can not take choice (iii), which is commonly used in the BR localization, in the form of 
Eq. ({TBI) . For this reason in addition to the reason that the r s -ambiguity itself is small, we 
take choice (i) in the present study. 

As for the ground state wave function of target nucleus, it is desirable to be as realistic 
as possible, and to be calculated theoretically because we are planning to apply the formu- 
lation to unstable nuclei where no experimental data are expected. Therefore, we employ 
the Hartree-Fock (HF) calculation with the finite-range Gogny force [ijj as an effective in- 



teraction. In particular, the D1S parameter-set [16| is adopted, which is applied widely and 
successfully to many nuclear structure problems (see, e.g., Ref.|17||). The standard method 
to solve the HF equation with the Gogny force is to expand the single-nucleon wave func- 
tions in terms of the harmonic oscillator basis. It is, however, not very accurate when the 
wave functions extend far outside nucleus due to the weak-bin ding , which is characteristic 



in unstable nuclei. The Gaussian expansion method (GEM) 
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19( is a powerful method to 



treat such a problem of the spacially extended wave functions, and it has been applied for 
solving the Hf[2(3] and HFB (Hartree-Fock-Bogoliubov) 21] equations. We have developed 



our own program to solve the HF and HFB equations by the Gaussian expansion, where the 
merit of the Gaussian form of the interaction is fully utilized. 

In the present work, the target nucleus is 90 Zr, which is a stable nucleus and the method 
of the harmonic oscillator basis expansion works. However, the use of the Gaussian basis 
expansion is still preferable because the calculation of the matrix elements of the ^-matrix in 
Eq. ({TBI can be done easily and accurately. We have calculated the ground state of 90 Zr by 
the HFB method. It is found that the neutron pairing gap vanishes because of the N = 50 
shell closure and the proton pairing gap is also very small due to the subshell Z = 40. The 
energy gain by the pairing correlation is less than 100 keV, and its effect on the density p v {r) 
is less than 1.5%. Therefore, we use the HF wave function neglecting the pairing correlation 
in this work. In more detail, the radial part of the single- nucleon wave functions in each 
iv] Zj)-channel are expanded by the Gaussian functions, 

n g 

0^) = E<? J)e ~ (r/Al)2 ' ( 14 ) 

1=1 

where we take n g = 14 and their ranges A ? ; (i = 1, .., n g ) are chosen to be from 1 to 5 fm by 
geometric progression according to GEM 18|, [l9| ; they are an almost optimal choice in the 
case of n g = 14. The coefficients C^f are determined by the HF variational equation. The 
resultant binding energy of 90 Zr is 785.995 MeV, which is compared to the experimental value 
783.894 MeV. This result corresponds to the calculation employing 25 shells (N osc < 24) in 
the harmonic oscillator basis. We believe that the obtained wave function is realistic enough 
to perform the test of the localization of optical potential. 



III. THE BRIEVA-ROOK LOCALIZATION 

A local potential [/^(R) trivially equivalent to the nonlocal potential U EX (R, r) is de- 
fined by 

U*f(R) X K,n (R) = / £/ EX (R r)xK,, (r) dr. (15) 

Brieva and Rook derived an approximate form (R) to the equivalent local potential 
U^(R) [1]. The derivation is composed of three approximations. The first approximation, 
called the local semi-classical approximation (LSCA) [22|, is 

Xk,ki W = Xk,ki (R + s ) ~ jrf££ (r) = xk* (R) e iK(R) - s , (16) 
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where the local momentum KK (R) is parallel to the flux of the scattering wave at R and 
its magnitude is determined to satisfy 

KK(R) = [2fi(E - V c (R)5»_\ /2 - t/ DR (R) - U^(R))]^, (17) 

i.e., KK (R) is evaluated self-consistently. LSCA has been successfully applied to studies on 
cross sections and spin observables for multistep direct (p, p'x) and (p, nx) processes as well 



as hyperon production cross sections 23]. 
LSCA yields a local form of 

^lsc(R) = E / <°^ R > r ^(s; pM)^ s ds. (is) 

The local potential C^lsc(R) of Eq. (|18|) is a function of the radial component R and the 
angle 6 between vectors R and K(R), as shown later. Obviously, LSCA is getting better as 
E increases. Actually, LSCA is good for E > 65 MeV, as shown later in Sec. HVl 



The second approximation, called the local Fermi-gas approximation (LFGA) 24j, is an 
approximation to the single-particle mixed density p„ 2 (R, r): 

p„ 2 (R, r) « p V2 (r m ) 3 [sin« 2 (r 3 )) - skl 2 {r g ) cos« 2 (r 9 ))] = p^ G (R, r), (19) 
where k^ 2 (r g ) is related to p V2 {r g ) as 

(F * 



/^ = %5-- (20) 
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LFGA is known to be a good approximation for small values of s 24j. The third approxi- 
mation is expressed by 

^ R > s ^ 3o (K(R)s), (21) 

where jx is the spherical Bessel function. This approximation is good when jo(K(R)s) <C 
jx(K{R)s) for X > 1. This condition is well satisfied when K(R)s < 1. In Eq. (TT8T) . 
the range of the integrand is about 0.5 fm because of the presence of the short-ranged 
interaction g^. In the surface region of A that is important for forward NA scattering, 
K(R) approximately equals the asymptotic wave number K. Hence, Eq. (l2Tj) is good at 
least for K < 2 fm -1 (E < 80 MeV). Eventually, the BR- type equivalent local potential 
£/br (R) is obtained by 

UE£(R) = £ / ^ FG (R> r)9^(s; Pv 2 (r 9 ))j (K(R)s)ds. (22) 



In the above derivation, LSCA is good for E > 65 MeV, while Eq. (12T|) is applicable 
for E < 80 MeV. However, as shown in Sec. IIV Al the BR-type local potential Ubr(R) = 
U BR (R) + U^(R) gives almost the same elastic-scattering cross section as that obtained 
by the exact calculation for 65 < E < 800 MeV. This means that the above derivation is 
not sufficient. Actually, as shown in Sec. IIV Eq. (I2"2"j) is derivable from Eq. (j!8p without 
the approximation (12TT) . that is, the approximation is redundant. Although, in the original 
work 1] of Brieva and Rook, the local momentum KK (R) is assumed to be parallel to the 
asymptotic momentum KK, it is not necessary because Eq. (122]) does not depend on the 
direction of K (R). Since LSCA itself is more accurate for KK (R) parallel to the flux of the 
scattering wave at R, we should think that the direction is also taken in the BR localization. 



IV. RESULTS 



A. Proton elastic-scattering from Zr 

Figure[Hpresents the differential cross sections of the proton elastic scattering from 90 Zr at 
(a) E = 65 MeV, (b) 185 MeV, (c) 400 MeV and (d) 800 MeV. For each panel, the horizontal 
lower (upper) scale represents the transferred wave number q (the scattering angle 6 cm ). The 
solid curves represent results of the exact calculation in which Eq. (CQ) is solved numerically. 
In the dashed curves, the medium effect is switched off from the exact calculation by taking 
g (s; p V2 = 0) in Eqs. (T5|) and ([3]). In the dotted curves, the exchange effect is neglected 
from the exact calculation by setting — in Eq. ([I]). Thus, the exchange effect is large 
at least up to E = 800 MeV, and the medium effect is significant up to E = 400 MeV. 

Now, the validity of the BR localization is tested. Figure [2] presents the same quantities 
as in Fig. [Tj but with different calculations. The solid curves represent results of the exact 
calculation of Eq. (CQ), while the dashed curves do results of the BR-type local potential 
Ubr(R), that is, the Schrodinger equation with Ubr(R) = U BR (R) + U^(R) is solved 
numerically. Seeing the difference between the two types of lines around q = 3.5 fm _1 , 
one can find that the error of the BR localization is getting small as E increases. For 
E = 65 MeV, the error is small at q < 1.7 fm _1 where the data are available, although it is 
sizable at large q around 3.5 fm~ . Thus, the BR localization is good for 65 MeV < E < 
800 MeV. We discuss this point in Sec. IIV Bl in detail. 
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In Fig. [3j the exact wave function (solid curve) of the Schrodinger equation Eq. ([I]) is 
compared with the approximate wave function calculated with Ubr{R) for the case of the 
p+ 90 Zr scattering at L = 8, where the elastic partial cross section af = 7r(2L+l)|l — Sl\ 2 /K 2 
calculated with the elastic S'-matrix element Sl becomes maximum. For both E = 65 MeV 
and 185 MeV, the approximation wave functions are very close to the exact ones. 

B. Validity of the BR localization 

First, we show the numerical test of LSCA, together with a simplified version of LSCA 
in which the absolute value of the local momentum, hK(R), is replaced by that of the 
asymptotic momentum UK. This version is referred to as LSCA-A in this paper. LSCA-A 
is, if justified, very useful since LSCA-A makes it much simpler the numerical task to obtain 

For this purpose, we consider the potential scattering and take Ubr(R) = U DK (R) + 
Ubr{R) as the potential. The exact wave function x~k,ui (r) of the potential scattering is 
compared with the approximate wave functions xWu ± ( r ) an d xWu± A ( r ) based on LSCA and 
LSCA-A, respectively. The wave functions are invariant under the rotation around the z 
axis, and hence the azimuthal angle ip of vector R can be set to zero. As an example, vector 
R is fixed to (R, 8, p) = (5 fm, ir/3, 0), and vector s in Eq. (fT5"|) is varied in a direction either 
parallel or perpendicular to R. For convenience, a variation of vector s in the direction 
parallel (perpendicular) to R is denoted by s r (sg); precisely, s r = (s • n)n with n = H/R 
and Sg = s — s r n. Note that LSCA is expected to work well if the potential around R varies 
slowly within the wave length of Xk,^i (*")• Thus, the choice of R = 5 fm can severely test 
the validity of LSCA. 

Figure H] represents the exact and approximate wave functions of the p+ 90 Zr elastic 
scattering at E — 400 MeV. In panels (a) and (b) where sg is varied with s r fixed to 
0, the wave functions Xk^ (r) (dashed curves) and xk S S" A ( r ) (dotted curves) agree with 
Xk,vi ( r ) (solid curves) within the range of the g-matrix interaction, i.e. at s < 1.5 fm. 
This is the case also for panels (c) and (d) where s r is varied with sg fixed to 0. Further, 
Xk S £ (r) is identical to A ( r ) at s < 1.5 fm. Thus, LSCA and LSCA-A are accurate 
for high E. At slightly lower energy of 185 MeV, as shown in Fig. \5\ the accuracy of LSCA 
and LSCA-A is almost the same as at 400 MeV. Figure O shows the results at 65 MeV. 
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One sees that LSCA and LSCA-A still work well, although its accuracy becomes slightly 
worse than for E > 185 MeV. Thus, we conclude that LSCA and LSCA-A are applicable for 
E > 65 MeV. Actually, as shown in Fig. [7] the elastic-scattering cross section for p+ 90 Zr at 
65 and 185 MeV calculated with LSCA-A in the BR localization agrees well with the result 
of the standard BR calculation with LSCA; particularly at 185 MeV, the difference between 
the two is within the thickness of lines. 

Next, we consider the remaining two approximations, i.e., LFGA and Eq. ff2T]) . One can 
find, however, the latter is redundant, if an extended LFGA below is justified. We consider 
the following approximation for the mixed density: 

p V2 (R, r) « p U2 (R) {skF \ R))3 [Mskl 2 (R)) - <(J2) ooB«(fl))] = pJf^R, r). (23) 

This approximation is referred to as LFGA-i?. Note that p^J G ~ R (R,r) is nothing but 
pJ;J G (R, r) of Eq. ffT9l) . with r m and r g replaced by R. Obviously, as long as LFGA is 
valid, LFGA-i? is also good for R larger than the range of the NN interaction s ~ 0.5 fm. 
Inserting Eq. ([23]) into Eq. flTBJ leads to Eq. (|22J without Eq. §23$, since p^ 2 (R,r) is a 
function of R and s in LFGA-i?. Thus, Eq. (I2lj) is not necessary, when LFGA-i? is taken. 

The validity of LFGA-i? is evaluated by comparing the results of the elastic-scattering 
cross sections calculated with LFGA and LFGA-.R. It is found that in the energy region 
of 65 MeV < E < 800 MeV, the two calculations show a perfect agreement. Thus, we 
conclude that Eq. ()2ip . which has imposed an upper limit of E where the BR localization is 
accurate, is actually redundant because of the good accuracy of LFGA-i?. This comes from 
the fact that, as mentioned above, LFGA itself is a very good approximation to the mixed 
density. In fact, it turns out that the p+ 90 Zr elastic-scattering cross section calculated with 
LFGA agrees very well with that obtained with explicitly using the mixed density. Another 
point to be mentioned here is that the simplest formula, Eq. fl20l) . is used for evaluating k^ 2 
from the one-body density. Our finding clearly shows that Eq. (1201) is enough to study the 



elastic-scattering cross sections; the higher-order corrections to k^ 2 [24], |29( are not necessary 
for this purpose. 

Finally, we comment on the r 9 -ambiguity of the BR-type local potential. Figure M shows 
the elastic-scattering cross section from 90 Zr at (a) E = 65 MeV, (b) 185 MeV, (c) E = 
400 MeV and (d) 800 MeV. The solid lines represent results of the BR local potential 
calculated with choice (i), and the dashed and dotted lines correspond to results of choices 
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(ii) and (iii), respectively. The differences among the three cases are appreciable only at 
dips for E = 65 MeV. For E = 185 MeV, they are appreciable even for tops. For higher E 
such as E > 400 MeV, the differences become negligible, since so is the medium effect itself. 
Thus, the folding potential has an appreciable r 9 -ambiguity only around E = 200 MeV in 
which the imaginary part of Ubr{R) is rather weak compared with the case of other E. 

V. SUMMARY 

We test the BR localization of the microscopic nucleon-nucleus optical potential over the 
wide range of 65 < E < 800 MeV and conclude that the localization is valid there. The BR 
localization is composed of the local semi-classical approximation (LSCA), the local Fermi- 
gas approximation (LFGA) and Eq. (j2ip . but these approximations can be reduced to two, 
LSCA and LFGA-i? (a modified version of LFGA). The former is reliable at E > 65 MeV, 
while the latter is good for any E. The approximate wave functions calculated with the 
BR-type local potential are very close to the exact ones. Thus, the BR-type local potential 
is quite useful in many applications, for example, as potentials between A and constituents 
of weakly bound or unstable projectiles. 
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FIG. 1: (color online) The differential cross sections of the proton scattering from 90 Zr at (a) 
E = 65 MeV, (b) 185 MeV, (c) 400 MeV and (d) 800 MeV. The solid curves represent the results 
of the exact calculation, while the dashed (dotted) curves corresponding to the calculation without 
the medium (exchange) effect. The horizontal lower (upper) scale shows the transferred wave 



number q (the scattering angle cm ). Experimental data are taken from Refs. 
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271, 



28]. 
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FIG. 2: (color online) Same as in Fig. [T] except that the dashed lines show the results of the BR 
folding model, i.e., with Ubr(R)- 
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the case of L = 8. The left and right panels represent the real and imaginary parts of the wave 
functions, respectively. The upper (lower) panels correspond to E = 65 MeV (185 MeV). 
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FIG. 4: The exact and approximate wave functions of the p+ 90 Zr elastic scattering at E = 
400 MeV. Vector R is fixed at R = (R,9) = (5 fm, vr/3). In panels (a) and (b) sq is varied, 
while in panels (c) and (d) s r is varied. The left (right) panels represent the real (imaginary) parts 
of the wave functions. In LSCA, the direction of vector K(R) is assumed to be pararell to the flux 
of the scattering wave at R. The same assumption is made also for the direction of K in LSCA-A. 
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FIG. 5: Same as in Fig. H but at E = 185 MeV. 
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FIG. 6: Same as in Fig. [1 but at E = 65 MeV. 
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FIG. 7: (color online) The validity of LSCA-A. Panels (a) and (b) represent the differential cross 
sections of the proton scattering from 90 Zr at E = 65 MeV and 185 MeV, respectively. The dashed 
curves denote the results of the BR local potential in which LSCA-A is taken instead of LSCA, 
while the solid curves show the results of the ordinary BR local potential. 
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FIG. 8: (color online) The r 9 -ambiguity of the BR local potential for the elastic-scattering cross 
section from 90 Zr at (a) E = 65 MeV, (b) 185 MeV, (c) E = 400 MeV and (d) 800 MeV. The solid, 
dashed and dotted curves stand for the results of the BR local potential calculated with choices 
(i), (ii) and (iii) for r g , respectively. See the text for details. 
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